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DETERMINATION OF GAS TEMPERATURES FROM 
LASER-RAMAN SCAHERING 


by Jack A. Saizman, William J. Masica, 
and Thom A. Coney 

Lewis Research Center 
SUMMARY 

An investigation was conducted on the feasibility of using Raman scattered laser 
light to measure local gas temperatures. Temperature errors are predicted based on 
the statistics of the scattering process and Raman spectral theory. Methods were de- 
veloped for the calculation of temperature measurement accuracies as a function of the 
Raman signal for different spectral line ratio schemes. 

Rotational Raman data were obtained from a laboratory- constructed Raman spec- 
trophotometer to provide experimental temperature measurements of varying accu- 
racies. Nitrogen was used as a gas sample for a range of controlled temperatures from 
253 to 313 K and pressures from 0. 5 to 5. 0 atmospheres. The temperature results 
from this Raman line intensity data indicate that the relations developed in this investi- 
gation can be used to predict temperature measurement accuracies. If the temperature 
accuracies indicated by these relations are adequate, real-time temperature measure- 
ments using Raman scattered light are feasible. Both the theoretical and experimental 
results of the study demonstrated that temperature measurement accuracies of ±7 K are 
possible within the limits of this investigation. 


INTRODUCTION 

Among the potential applications of the vmique properties of laser light are those 
techniques that utilize the laser as a remote sensii^ probe. A number of these techni- 
ques exploit the laser as an optical radar device (lidar) whose return signal results 
from the backscattering of the laser light by gas molecules. (A comprehensive survey 
of these techniques is contained in ref. 1. ) Investigators have used the elastically back- 
scattered light from a laser beam to remotely observe a number of meteorological 



variables such as aerosol or particulate concentration levels (refs. 2 and 3). An inher- 
ent limitation of this probing method results from the fact that the return signal is the 
combined result of Rayleigh scattering by molecules and Mie scattering by aerosols. 

And there is no convenient method available to separate these two effects since they both 
return at the transmitted light wavelength. Also, because all the received information 
is contained in intensity variations of the transmitted wavelength return, extraneous ef- 
fects of the intervening atmosphere on this return cannot be determined exactly and must 
be calculated by using approximations. 

More recently, it has been suggested that lidar techniques and even some short- 
rai^e devices could utilize Raman scattered light to provide additional data on meteoro- 
logical state parameters. The Raman return signal is the result of inelastic scatterii^ 
processes which occur when the photons of the incident light beam interact with the gas 
molecules and the energies of a small portion of the scattered photons are either in- 
creased or decreased relative to the exciting photons. This chaise in energy and the 
resulting wavelength shift in the scattered light occurs in quantized increments which 
correspond to differences in the energy levels of the scattering molecules. The Raman 
spectral line positions are independent of the exciting wavelength and scattering angle 
and are dependent only on molecular structure. The Raman return signal can, in theory, 
yield information on atmospheric pressures, temperatures, and molecular species con- 
centrations. The Raman technique has already been employed to differentiate the lidar 
signal return from a number of basic atmospheric constituents and possible pollutants 
(refs. 4 to 6). 

Utilization of the Raman return of a lidar system to remotely measure local atmo- 
spheric temperatures was initially proposed by Cooney (ref. 7) and expanded by Press- 
man (ref. 8). The method involves simultaneously monitoring at least two portions of 
the Raman return signal from a major atmospheric constituent and examining the ratios 
of their intensities. An especially attractive feature of this technique is that, if the fre- 
quencies of the measured Raman signals are relatively close, the effects of all the ex- 
traneous variables (such as intervening atmospheric phenomena) are essentially elimi- 

Q 

nated. Unfortunately, the reduced intensity of the Raman signal (i. e. , at least 10 lower 
than the Rayleigh signal) has hindered actual lidar temperature probing. In fact, no 
Raman-based temperature data are available for near-ambient atmospheric gases. 

This report presents the results of an investigation on the feasibility of using Raman 
scattered light to measure local gas temperatures. Temperature errors are predicted 
based on the statistics of the scattering process. Methods are presented for the calcu- 
lation of temperature measurement accuracies as a function of Raman signal return 
magnitude for different spectral line ratio schemes. These theoretical calculations are 
compared with experimental data obtained from a laboratory- constructed Raman spec- 
trophotometer. Nitrogen was used as a gas sample for a range of controlled tempera- 
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tures from 253 to 313 K and pressures from 0. 5 to 5.0 atmospheres. Only pure rota- 
tional Raman spectral lines were considered. 


SUMMARY OF RAMAN THEORY 

The theory of Raman scattering and its application to molecular analysis is well 
known. Reference texts on the Raman effect include those by Bhagavantam (ref. 9), 
Hertzberg (ref. 10), Stoicheff (ref. 11), and Szymanski (ref. 12). This section briefly 
reviews pertinent Raman theory and the application to remote sensing of thermodynamic 
variables. 


Raman Theory 

The Raman effect is a quantum mechanical process wherein incident light quanta ex- 
change energy with a molecular scattering site. Exchanges can occur with rotational, 
vibrational, etc. , molecular energy levels, resulting in a spectrum of scattered light 
ordered by the energy levels in the molecule. The intensity per molecule of one Raman 
line associated with an energy transition between an initial molecular state m and a 
final state n is given by 


where is the wavenumber of the incident radiation and is far removed from an ab- 

° 2 
sorption band, is the wavenumber of the scattered radiation and is the 

probability of transition. The wavenumber shift Ap is defined by the transition energy 

difference 






m 


he 


( 2 ) 


and can be greater than, less than, or equal to zero depending on the relative magnitude 

of the molecular energy levels. The Raman line for > 0 is called the anti-Stokes 

line; the reverse transition, Ai> < 0, is called the Stokes line. If there is no net energy 

exchange during the scattering process. Ay = 0 and the scattered radiation is the fa- 

2 

miliar Rayleigh scattering. The transition probability is obtained from consider- 

ations of the dipole moment induced in the scattering molecule by the incident radiation. 
In general, these probabilities are functions of the symmetry properties of the scattering 
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molecule; the energy states involved in the transition; and the intensity, direction (rel- 
ative to the angle of observation), and polarization of the incident radiation (ref. 12). 
For a given state of incident radiation and a given scattering system, is a function 

solely of molecular structure. 

If N is the number of molecules in the initial state, it is assiuned that the total 
m 

intensity of a Raman line is given by 


= N^i _ 
mn m mn 


(3) 


The probability of occupation of a state with energy 
tical thermodynamic expressions of the form (refs. 
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-E/kT 
e 

All states 


can be calculated from statis- 
m 

13 and 14) 

-E^/kT 

= £ (4) 

q 


Substituting, 


Ne 


-E/kT 
m' 


ran 


mn 


(5) 


The quantity q is called a partition function and is ordered to the approximated molecu- 
lar energy level structure; for example, there are rotational partitional functions de- 
fined by rotational energy levels, etc. (ref. 13). Since i^^^ is characteristic of the 
molecular structure of the scattering site, a Raman spectrvim measurement can, in 
principle, provide composition identification. Once the composition is specified, meas- 
urement of could yield temperature and density data. These considerations form 

the basis for using Raman scattering as a remote sensor. 


Temperature Data From Raman Scattering 

A potential application of remote Raman sensing is obtaining temperature data in 
the troposphere. It is known that the pure rotational nitrogen lines are the most intense 
lines in a Raman spectrum of air near normal temperature and pressure (NTP) condi- 
tions and, therefore, these lines were selected for consideration in this study. 

The intensity of a pure rotational Raman line for a linear molecule can be expressed 
as 
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= CqNj(i;q + Ai/)‘^| (6) 

where is an intensity factor that is a function only of rotational quantum numbers and 
Cq contains those parameters that are constants for a fixed scattering system. The 
term Cq is a function of the index of refraction at the incident and scattered wavenum- 
bers (ref. 15) which includes, in the case of remote sensii^, atmospheric effects on the 
incident and return radiation. For pure rotational Raman scattering in gases, the dif- 
ference between the incident and scattered wavenumbers is less than 200 cm“^ (less than 
1 percent in the visible optical spectrum). Therefore, the indices of refraction are ef- 
fectively equal and Cq is independent of Ay. 

In rotational Raman scatterii^, only certain transitions between rotational energy 
levels are allowed. The selection rules for a linear molecule are j - Z =0, ±2, and the 
corresponding Stokes and anti-Stokes intensity factors are 


b| 


j+2 


3(j + l)(j + 2) 
2(2j + l)(2j + 3) 


b? 3j(j - 1) 

^-2 2(2j + l)(2j - 1) 


(7) 


The rotational energy levels of diatomic molecules can be adequately approximated 


by 


E. = j(i + l)k0 


( 8 ) 


where 0 is a constant, related to the moment of inertia of the molecule, called the ro- 
tational characteristic temperature. From eqviations (2) and (8), the Raman displace- 
ments are 



The constant 0 can be determined from spectroscopic measurements of Ay and has a 
value for nitrogen of 2. 86 K (ref. 13). 
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The rotational partition function for a diatomic homonuclear molecule with even 
mass number is 


^rot 


_ g(g + 1) 


j=0,2,4, 


(2j + 


2 j=l,3,5, . . . 

where g = 2 t + 1, Th/2n is the nuclear spin, t = 1 for nitrogen. Then 


^ . g(g ± l)(2j l)e 

N " 2Vot~ 

where +1 is chosen for even j, -1 for odd j. Substituting equations (7), (9), and (10) 
into equation (6) yields the intensity of a Stokes Raman line 




CjNg(g ± l)(j + l)(j + 2) 


L 

j 

[ ° he \ 


j + ^Kot 



,-j(j+l)0/T 


and an anti-Stokes line 


-AS 
■] 


CjNg(g ± l)(j)(j - l)j^i^Q |j 

(2j - Dq^ot 




( 11 ) 


( 12 ) 


where Cj is a new scattering system constant assumed to be independent of j, as pre- 
viously noted. The development has assumed a diatomic homonuclear gas structure 
such as nitrogen, but expressions similar in form are obtained for other molecular 
species. The transition j — j - 2 in equation (12) demands that j > 2 and that the first 
allowed anti- Stokes transition be from j=2 to j = 0. Conventionally, the first Raman 
line on either side of the exciting line is labeled 0, corresponding to the lower energy 
level in the transition. If J = 0, 1, 2, . . . denote Raman line position, then J = j for 
Stokes' lines and J = j - 2 for anti-Stokes' lines. Equations (11) and (12) become for 
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nitrogen 


^ CgNCJ + 1)(J + - 7. 96^J + 


,-J(J+l)0/T 


(13) 


rAS 


C2 N(J + 1)(J + 2) I^i/q + 7. 96^J + e 


,-(J+2)(J+3)0/T 


(2J + 3)q^^t 


(14) 


and I, 


where C 2 contains the factor g(g ±1). If the intensities of any two lines Ij 

are ratioed, all the factors dependent on the scattering system are canceled, including 
atmospheric path effects in remote sensing, leaving 


h. -K(Ji,J2)0/T 

= C(Jp J2)e 


or 


T = 


K(JpJ2)0 

7cWi.J2»J2, 


(15) 


where K(Jj^, J 2 ) and C(Jj^, J 2 ) are known from equations (13) and (14). Therefore, it is 
theoretically possible to calculate temperature from the remote measurement of two 
scattered Raman rotational line intensities. 


RESULTS AND DISCUSSION 
Temperature Measurement Error Calculations 


In practice, the feasibility and accuracy of the proposed method of temperature 
measurements depend on how the Raman line intensity ratio varies with temperature 
and how accurately the ratio can be measured. If a(Ij /ij ) is the error in the 
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measurement of Ij and S Is the variation of the line intensity ratio with respect 

to temperature, the error in the measurement of the temperature is 



The variation of the Raman line intensity ratio (eq. (15)) with respect to temperature 
is 


S = 




Maximizing S by the proper choice of Jj and 3 ^ yields the most sensitive measure- 
ment over a given temperature range. For 0/T values of the order of 10 S mono- 
tonically approaches a maximum with increasir^ values of Jg for values of Jj around 
4 to 6. Examples are shown in figure 1 for Jg equal to 16 and 20. Unfortunately, the 
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(a) Stokes -to -Stokes ratio. (b) Anti -Stokes to anti -Stokes ratio. 

Figure 1. - Variation of line Intensity ratios. 
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Raman line intensities rapidly decrease at large values of J, where the exponential term 
in equations (13) and (14) dominates. Increasing S ultimately causes the accuracy of 
the temperature measurement to decrease because of the difficulties in accurately meas- 
uring the intensities of the higher J-value Raman lines. The optimum selection of the 
Raman line is, therefore, primarily determined by an analysis of the line measurement 
error. 


Optimum selection of I 


■i/’ 


ties directly affect the error A 
perature calculation. The total error A 


It . - Measurement errors in the individual line intensi- 

J2 

involved in the intensity ratio used in the tern- 






results from the combination of exper- 


imental error and the statistics of the scattering process. It is possible to estimate an 
optimum temperature accuracy (i. e. , minimum AT) for the inherent scattering error. 
When considering measurements of Ij in terms of events per unit time (such as emitted 
photons or recorded counts per second), this inherent error can be estimated by a stand- 
ard Poisson distribution approximation. Since the recorded counts per second or ob- 
served count rate includes both a Raman signal count Ij and a bacl^round count rate Ig, 
the standard deviation of Ij is (Ij + 2Ig)^/^. Then, the standard deviation of the inten- 
sity ratio is calculated to be 


a 


r,B - 



and if Ig « Ij, 


a 


r 



(18) 


(19) 


From equations (19) and (16), a minimum standard deviation of the calculated tempera- 
ture can be written as 


T ■ 



0K(Jj,J2) 



1/2 
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and, if the root-mean- square error is accepted as the best error representation 


= AT = 


0K(Jj, Jg) ' 


^4 ;^2 

dj jg 


1/2 


( 20 ) 


In this report, arp and AT will be used synonymously. If another error representation 
were chosen, AT would be equal to multiplied by a constant factor. 

Some a^p's at 300 K are listed in table I to illustrate their dependence on Jj, Jg, and 
signal coimt rate. These values were calculated from equation (20) by assuming a coimt 
rate for the J =6 Stokes Raman line and then using equation (15) to obtain the count rate 
for other lines. The Stokes J =6 line is the most intense Raman line in nitrogen near 
NTP conditions, and for this reason was selected for comparison of signal magnitudes. 

O 

The count rates for Ij_g in table I approximate those values recorded experimentally 
at different gas pressures. 

g 

With the Raman signals expressed relative to lj_ 0 , O'rp can be written in the form 


CT 


T ~ 


M(J, T) 



( 21 ) 


and the least error in temperature measurement corresponds to minimum values of the 
function M and, obviously, the largest absolute signal strength possible. These mini- 
mum values of M at three temperatures are given in table II. These values yield the 
optimum temperature measurement accuracies that are possible. The corresponding 
values of and J 2 are, therefore, optimum Raman line numbers. The error pre- 
dictions presented in tables I and II indicate that the minimum temperature measurement 
error will vary from ±4 to ±15 percent with Raman signals raiding from 500 to 50 coimts 
per second. 

Multiple line ratios . - In an attempt to enhance the temperature measurement accur- 
acies previously obtained, a calculation technique was examined that employed ratios of 
the summed intensities of a number of Raman lines. The smnmation was restricted to 
adjacent spectral lines to approximate a spectral wavelength interval (e. g. , as might be 
seen by an optical bandpass filter). The ratio of two sums of Raman lines is 
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( 22 ) 


P = 


z 


1 



2 


where Ij has the form (eqs. (13) and (14)) 


Ij = C3f(J)e 


-H;j(0/T) 


One technique of solving for the temperature from equation (22) is to expand the expo- 
nential to the first order and use this as the basis for numerical iteration. For the tem- 
perature range 253 to 313 K, 0/T varies from 0. 01130 to 0. 00913 for nitrc^en. Let 

0 

— = a + X 
T 


such that 


-&t(®/T) 

e = e e (1 - ^jx) 


(23) 


which is accurate to better than 12 percent with J < 20 and a = 0. 010215. Substituting 
into equation (22) 


^f(Jl)e (l - 


^f(J2>e ^■’2 (l - 








(24) 


^f(J2>e ■’2 
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Note that, as with the previous ratio of sir^le Raman lines, all the terms dependent on 
the scattering system are eliminated. Solving for T, 



Using a value of T from equation (25), one can expand the exponential in equation (23) 
in a Taylor series about T to improve the calculation. Over the range of variables in 
this study, only three iterations (at most) are necessary for satisfactory convergence. 

An analysis of the expected temperature measurement accuracy for the technique of 
using ratios of summed line intensities was performed similar to that previously dis- 
cussed for ratios of single line intensities. The ratio response S at discrete tempera- 
tures was determined by computing the value of the ratio p at 2. 5 K on either side of 
the specified temperature and calculating the slope of Ap against AT. In estimating 
the ratio measurement error a^, it was assumed that the errors of the individual line 
intensity measurements were independent and, hence, additive in the summation proc- 
ess. Selected results of the calculation of Orri from a and S for a Raman Stokes 

O J- P 

spectrum are presented in table HI for an Ij_g = 500-counts-per-second signal level. 

In order to limit the number of possible ratios and facilitate the calculations, only 
even-numbered J lines were considered in the summations. If the odd-numbered lines 
are included between the adjacent even lines, the ratio measurement error is de- 
creased, while the slope remains nearly unchanged. These effects result in a decrease 
in Orp, the temperature measurement error, when the odd lines are included. In fact, 
the inclusion of the lower-intensity odd-numbered J lines within the spectral interval 
will reduce the temperature measurement error by as much as 10 percent from those 
values of listed in table HI. 

Comparing the results in table HI with those in table I indicates the desired improve- 
ment in temperature accuracy if spectral intervals, rather than discrete lines, are used 
in the calculation schemes. The Orp of 6. 9 K as presented in table III is the minimum 
at that signal return and temperature for all the summation ratio schemes (i. e. , Stokes 
to Stokes, anti-Stokes to anti-Stokes, and Stokes to anti-Stokes). Over a return signal 
range from 500 to 50 counts per second, the minimimi temperature measurement error 
range at 300 K is from ±2. 3 to ±7 percent for this summation technique. The particular 
width and position of the spectral regions which give an optimum temperature accuracy 
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at other temperatures can be obtained from calculations similar to those used in tabu- 
latii^ table HI. 


Experimental Results 

To establish the applicability of the preceding relations and to confirm the tempera- 
ture measurement error calculations, a series of experiments were conducted using a 
Raman laser spectrophotometer. The instrument and experimental technique are de- 
scribed in appendix B. 

Feasibility of Raman method . - The first set of experiments were designed to an- 
swer the basic question of whether in fact one could experimentally measure tempera- 
ture from Raman spectral lines at the temperatures and pressures of interest. A series 
of tests were run with the nitrogen sample at 1-atmosphere pressure, with three repeti- 
tions each at 253, 263, 273, 283, 293, 303, and 313 K. The order in which the tests were 
run was completely randomized. Data retrieved from these tests were obtained with an 
analog recording system. A portion of a typical spectrum is shown in figure 2. 

The intensity alternations between the odd and even J lines seen in figure 2 are as 
predicted by the statistical weights given in equation (10). Raman line intensities were 
obtained by subtracting an averaged background and noise count rate from each of the 



Figure 2 , - Analog data plot of rotational Raman spectrum of nitrogen. 
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analog peaks. For this test series, the average maximxim line intensity (i. e. , the 
Stokes J=6 line intensity) was of the order of 120 counts per second. 

All possible single-line to single-line intensity ratios formed with these data were 
subjected to a linear regression analysis. The model employed in the regression com- 
putations was of the form (cf, eq. (15)) 

= (y) - =5 

Each of the separate line ratio correlations which resulted from using this model was 
subjected to an analysis of variance and other standard statistical tests. (Ref. 16 con- 
tains a complete description of the tests. ) These tests were employed as selection 
criteria in a screening process used to reduce the more than 800 possible ratios to a 
smaller number where there was significant and sufficient correlation between the de- 
pendent and independent variables. The first selection test resulted from an analysis 
of variance of the overall significance of the model. Only those ratio correlations which 
exceeded a severe F test (i. e. , F > 40 for one and 19 degrees of freedom) were se- 
lected for further analysis since this test indicated a significant amount of variance re- 
moved by the correlation as compared to the residual variance remaining. Also, be- 
cause of the replications available, these correlations were tested with an analysis of 
variance of the lack of fit to be certain that a linear regression was sufficient to account 
for the variation in the intensity ratio. 

The group of selected line ratios satisfying these first statistical tests were further 
screened by rejecting all those where the standard error of estimate exceeded 10 per- 
cent, The 32 line ratio combinations which remained after this screening process were 
then considered to have their variation with temperature sufficiently correlated to allow 
comparison with the predicted theoretical results. These experimental correlations 
were compared with equation (15) by comparing the computed correlation constants with 
those predicted in the equation. The average difference between the correlations and the 
predicted constants was less than 4 percent. 

The excellent agreement of these results show that the relation in equation (15) can 
be used to calculate temperatures from experimentally measured rotational Raman line 
intensities and their ratios. Figure 3 shows the results for one such ratio, the Stokes 
J=4 line to the Stokes J=16 line. Although, not predicted by theory, this ratio yielded 
the best temperature measurement accuracy for this series of tests. The theoretical 
results as listed in table n indicate that the ratios S2/S16 and S2/S18 should have 
given better temperature measurement results. It was found, however, that a large 
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Figure 3. -Temperature measurement results calculated from intensity ratio of 
Stokes J=4 line to Stokes J=16 line. 


first-order grating ghost was superimposed on the Stokes J=2 line and, thus, the in- 
tensity of this line could not be measured without a large degree of error. This situation 
was corrected for further testing by the use of a polarization rotation as discussed in 
appendix B. From this initial series of tests, then, it was not possible to establish con- 
clusively the optimum line ratio combination. 

Temperature me asurement errors . - A second series of tests were conducted to 
investigate the temperature measurement errors evidenced by the data spread in fig- 
ure 3. These tests were designed to examine how accurately the error predictions pre- 
viously developed represent actual deviations in experimental line intensity data and the 
corresponding calculated temperatures. Because these errors were predicted to be a 
function of the absolute Raman line intensities (i. e. , in events or counts per second), 
this series of tests was conducted such that the intensity of the scattered Raman light 
could be varied. The intensity was controlled by adjusting the gas pressure in the test 
cell and thus varying the molecular density in the scattering volume. All these tests 
were made at a gas temperature of 299 K; and four repetitions were made at each of the 
pressures 0. 5, 2. 0, and 5. 0 atmospheres, yieldir^ maximvim Raman line intensities 
(i. e. , Stokes J=6 lines) which average 37, 289, and 519 counts per second, respec- 
tively. 
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The analysis of these Raman line intensities was restricted to even-numbered lines 
in order to reduce the task of data acquisition and analysis and to limit the number of 
ratios in the investigation. It was felt that no loss in generality resulted from this re- 
striction because the ratio-temperature relation of the odd lines had been examined in 
the previous series of tests. The four repetitions at each of the different count rate 
levels were used to calculate a standard deviation for the individual line ratios. In 
evaluating the standard deviations, a calculated line ratio value (from eq. (15) evaluated 
at 299 K) was used as a mean value. Each of these standard deviations was then com- 
pared to a predicted error as calculated from equation (18) using mean line intensity 
values. 

In evaluating the analc^ data errors, it was found that the contribution due to back- 
ground error could be ignored, reducir^ equation (18) to the simpler form of equa- 
tion (19). Predicted errors calculated from equation (19) for the two extremes in count 
rate levels (i. e. , 37 and 519 counts per second) are compared with the measured errors 
in figure 4 for even Stokes line ratios. As shown in figure 4, the actual error involving 
all contributions was generally less than that calculated from equation (19), which con- 
siders only the inherent scattering distribution. This was because of the integrating 
properties of the resistance-capacitance (RC) circuit in the analc^ recording system, 
which resulted in the measurement of a count rate averaged over several seconds rather 
than a discrete 1- second value. 



Figure 4. - Analog intensity ratio data for Stokes even lines. 
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Figure 5. - Analog temperature data for Stokes even lines. 


Experimental temperatures were then calculated for this data using equation (15). 
Variance in these measured temperatures for each repetitious set were then compared 
with the calculated temperature errors predicted from equation (20). The results are 
presented in figure 5. As with the line ratio error, the temperature measurement error 
was consistently lower than that predicted. However, the experimental data did follow 
the functional trend of the predicted error and, hence, calculated values of ffrp provide 
at least conservative estimates of the temperature measurement accuracy. Although 
the results presented are only for Stokes ratios, they were found to be equally applicable 
to anti-Stokes-to-anti-Stokes ratios and Stokes-to-anti- Stokes ratios. 

A similar analysis was performed with the line intensity data in digital form. Dis- 
crete 1-second count rates for each Raman line were obtained by plottii^ the data as in 
figure 6 and measuring that count rate nearest the center of the peak contour. As shown 
in figure 7, the measured line ratio errors for the digital data were slightly higher than 
those predicted by using equation (18). This indicates the effect of the other error con- 
tributions, such as experimental and measurement errors, which accompany the vari- 
ance due to the statistics of the basic scattering process. Experimental temperatures 
were again calculated and the errors in these temperatures are compared in figure 8 
with those predicted. 


17 



CD 


60 r- 


50 


40 


30 


20 


10 


• Included in the line 
intensity calculation 


O 

( 5 ^ 

I — Q. 




-Contour center 




O f 

oT 


cP 

o 


o 

o 


o 

o 

o 




o 

o 


o 


99S^c^ 




o 

o 




Wave number, v, cm' 


Figure 6. - Raman line plot of digital data. 




Figure 7. - Digital intensity ratio data for Stokes even lines - 
counting for 1 second on each line. 




Figure 8. - Digital temperature data for even Stokes lines - 
counting for 1 second on each line. 


The actual errors (as shown In figs. 7 and 8) can be reduced to a value nearly the 
same or below those predicted in equation (18) by using a line count rate value which is 
averaged over several seconds. By averaging the peak count rate with the count rate 
values on both sides of it and using this value to calculate a ratio and a temperature, 
the measured error was less than that predicted for discrete 1-second coimting (figs. 

9 and 10). 

By averaging more count rates grouped around the center of the peak contour, the 
measured error can be reduced further (i. e. , by approximately the square root of the 
number of discrete 1-second count rates which are being averaged). In this investiga- 
tion, the spectrum scanning speed limited the number of 1-second counts that could be 
averaged to five. Those values which were more than 2 seconds from the contour cen- 
ter could no longer be considered peak count rate values but were contained in the wii^s 
of the contovir. 

Multiple-line ratios . - The same Raman line intensity data used for the single-line 
ratio calculations were also employed to examine the ratio-of-summed-lines calculation 
technique. The results of these calculations agreed favorably with the predictions. A 
comparison of the experimental temperature measurement error with the predicted 
error is shown in figure 11 for ratios of two, three, four, and five smnmed lines. The 
increase in accuracy predicted for the technique is illustrated in figure 11 by the num- 
ber of CTrp, both calculated and measured, that are less than 10 K. 
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Measured standard deviation of intensity ratio 



Calculated standard deviation of intensity ratio, g 


Figure 9. - Digital intensity ratio data for Stokes even lines - 
averaging over 3 seconds on each line. 


Measured standard deviation of temperature, 



Calculated standard deviation of temperature, a-p = Oj. g/S, 


Figure 10. - Digital temperature data for even Stokes lines 
counting for 3 seconds on each line. 


□□ 
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Calculated standard deviation of temperature, Oj = Op/S, K 

Figure 11. - Temperature data for ratios of summed line inten- 
sities. 


CONCLUDING REMARKS 

The purpose of this investigation was to examine the feasibility of usir^ Raman 
scattered light to measure local gas temperatures. The statistics of the scattering 
process provide a basis for error predictions. Methods developed using this error 
model and basic Raman theory can be used to calculate theoretical temperature meas- 
urement accuracies as a function of Raman signal return for different intensity ratio 
schemes. Sufficient experimental rotational Raman line intensity data were obtained 
with gaseous nitrogen at pressures from 0. 5 to 5. 0 atmospheres to provide temperature 
measurements of varying accuracies over a temperature range from 253 to 313 K. The 
results from both analog and digital line intensity data used in both single-line and 
summed-line ratios indicate that the relations developed can be used to effectively pre- 
dict temperature measurement accuracies. If the temperature accuracies indicated by 
these relations are adequate, real-time remote temperature measurements usii^ Raman 
scattered light are feasible. 

Although the summed line intensity data were obtained by scanning and summing dis- 
crete peak intensity measurements, the results imply that enhanced accuracies can be 
gained in real-time measurements if the total intensities of two spectral intervals are 
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used in the temperature calculations rather than the ratios of single lines. The particu- 
lar width and position of these spectral regions which yield an optimum temperature 
measurement accuracy can be obtained from calculations similar to those given in this 
report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 28, 1971, 

124-08. 


22 



APPENDIX A 


SYMBOLS 


cCJj.Jg) 

Cqj Cj, C 2 , Cg 


intensity factor 

constant determined by the particular line intensity ratios being con- 
sidered, e.g. , 

4 

scattering system constants, cm J sec" 
correlation constants 
speed of light, 3x10^® cm sec"^ 

energy of quantum states j, m, and n, respectively, J 



g 

h 

% 

^J’ ^n 

■! 

tAS 

^mn 

J 

j 

K(jj, Jg) 

k 

I 


spin degeneracy 

Planck’s constant, 6.625x10“ J sec 

bacl^round intensity, including stray light and electronic noise, 

_i 

counts sec 

total Raman intensity for the J and j lines and for the transition from 
level m to level n 

J.-U -I 

total Raman intensity of the Stokes line in J sec"^, photon counts per 
second 

total Raman intensity of the anti-Stokes line in J sec"'*^, photon 
counts per second 

rms error in measurement of line ratio 

Raman line intensity per molecule for the level m to level n transi- 
tion, J sec"^ 

Raman line designation, for Stokes spectral lines J = j, for anti- Stokes 
lines J = j - 2; J = 0, 1, 2, . . . 

rotational quantum number (initial energy state) 

constant determined by the particular line intensity ratios being con- 
sidered 

Boltzmann constant, k = 1. 381x10“^^ J K”^ 

rotational quantum number (final energy state) 
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M(J, T) 


m 

N 


n 


P 


2 

mn 


q 

S 


T 

AT 

0 


P 


Ap 


a 


r 


a 


r,B 


a 


T 


a 

P 

T 


function of J and T relating the standard deviation of the temperature 
measurement to the intensity of the J=6 Stokes line, K 

quantum nvimber (initial energy state) 

total number of scattering molecules 

quantum number (final energy state) 

transition probability, sec"^ 


partition function 

derivative of line intensity ratio with respect to temperature. 


S = 



-1 


temperature, K 

error in measurement of temperature, K 
rotational characteristic temperature, K 
wave number of exciting radiation, cm”^ 

wave number shift, from p^, of scattered radiation, cm 



error in ratio of two sums of Raman spectral lines, e. g. , 


Vj, 



standard deviation of ratio of lines 

standard deviation of ratio of lines including bacl^roimd effects 
standard deviation of temperature, K 


standard deviation of ratio of two sums of lines 


nuclear spin number 
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APPENDIX B 


APPARATUS AND EXPERIMENTAL TECHNIQUES 
Test Apparatus 

A laboratory Raman spectrophotometer was constructed to provide experimental 
data to compare with the theory. The instrument was used to make Raman spectrum 
measurements of gaseous nitrogen over a temperature range from 253 to 313 K at pres- 
sures ranging from 0. 5 to 5. 0 atmospheres. The spectrophotometer was built as a mod- 
ular system to increase instrument versatility. Basic units of the instrument are shown 
in the block diagram and accompanying photograph in figure 12. 

The exciting light source was a continuous-wave argon-ion laser which was tunable 
to a number of wavelengths through the use of an intracavity Littrow prism. For this 
investigation, it was operated at a wavelength of 4880 A, which provided approximately 
200 milliwatts of multimode power. The output beam was linearly polarized to better 

3 

than 1 part per 10 with the electric vector oriented vertically. 

An external optics arrangement (fig. 12) was employed to increase the scattered 
light radiance, as described in reference 17. After initially being directed toward the 
gas sample by a high- reflectivity mirror, in some tests, the laser beam was passed 
through a crystalline quartz plate which rotated the beam polarization by 90°. (This was 
done to minimize the Rayleigh light being scattered toward the collecting optics.) The 
beam was condensed by a focusing lens to create a high-power-density scattering site in 
the gas sample. The optical path of the beam was terminated at a return mirror which 
refocused the beam back into the scattering site, thus creating repetitive passes of the 
beam through the gas sample and increasing scattered light radiance. The volume of 
maximum light scattering or the "Raman light source" could be approximated by a hori- 
zontal cylinder 20 micrometers in diameter and 2 millimeters in length (ref. 17). In- 
cluded in the external optics was a light sensor which monitored the power of a secondary 
laser beam. The output of this sensor was observed on a photometer which provided 
continuous recording of the relative laser power. 

_3 

The gas samples were confined in a cell approximately 10 cubic meter in volume. 
Temperature regulation of the gas was provided by flowing liquid through the walls of the 
cell and using a controlled-temperature liquid bath. The temperature of both the cell and 
the gas, along with the gas pressure, were continuously monitored throughout each test. 
The optical windows in the cell through which the laser beam passed were parallel plates 
of coated quartz. At right angles to the axis of these windows were two optical ports 
which comprised a portion of the collecting optics system. At the lower gas tempera- 
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Beam optics 
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Detection and recording electronics 



Figure 12. - Raman spectrophotometer. 
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tures, all windows and ports were continuously purged with dry nitrc^en to prevent fog- 
ging. 

Light scattered at right aisles to the axis of the exciting laser beam was collected 
over a large solid angle by a lens which was a port in the cell. This lens, a field stop, 
and a secondary lens comprised an optical system (fig. 12) which magnified the image of 
the Raman source cylinder by a factor of 2. 5 and matched the optical acceptance angle 
of the monochromator. By rotating the magnified image through 90° with a prism, it 
was possible to efficiently fill the vertical slits of the monochromator with collected 
light. The intensity of the collected light was enhanced by gathering that light which was 
scattered 180° away from the primary collection lens with a concave mirror. Because 
of the polarization sensitivity of the monochromator, a quartz wedge depolarizer was 
placed immediately in front of the entrance slits of the monochromator. 

Wavelength discrimination was provided by a single 3/4-meter Czerny-Turner grat- 
ing monochromator. Both slits were operated at a nominal 2. 3 cm ^ (a mechanical slit 
width of 50 jim.) to provide both efficient coupling with the collecting optics and adequate 
wavelength resolution. For this investigation, the wavelength scanning speed was fixed 
at 2 A per minute. Because the monochromator was only a single-pass instrument, 
stray light was a problem if the intensity of the collected Rayleigh light and spurious re- 
flected laser light was very high. Fortunately, the Rayleigh light was highly polarized 
and its collected intensity was reduced by the 90° rotation of the exciting beam polariza- 
tion. With this procedure, the collected light at the excitation wavelength was suffi- 
ciently reduced to keep stray light and primary gratii^ ghosts at an insignificant level. 

Light emerging from the exit slit was focused down by a multielement lens onto a 
1-millimeter-diameter, S-20 surface of a photomultiplier tube which employed a chan- 
nel electron multiplier to provide gain. A highly stable power supply provided the 
multiplier bias, which was normally operated at 2250 volts. The small size of the photo- 
cathode surface and the low noise characteristics of the multiplier circuit made it pos- 
sible to operate the tube without cooling. When it was operated at room temperature, 
the combined noise and background count rate was of the order of 20 counts per second. 

The signal from the photomultiplier tube was sent through a system of photon 
counting (or pulse counting) electronics (fig. 13) for amplification and display. Pulse 
signals were first passed directly to a low-noise preamplifier for impedance matching 
and first -stage amplification before transmission to the remaining rack-mounted elec- 
tronics. A pulse-shaping linear amplifier then conditioned the signal for input into a 
pulse height analyzer. The pulse height analyzer served as a discriminator-trigger cir- 
cuit so that each signal pulse was weighted equally and converted into a narrow-width 
pulse. Coincidence outputs from the pulse height analyzer were channeled to two sepa- 
rate recording systems, one of which was an analog system and the other a digital 
system. 
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Figure 13. - Detection electronics. 


The analog system consisted of a variable RC network or ratemeter and a chart re- 
corder. The ratemeter was operated in the linear mode. The RC circuit was normally 
operated with a time constant of 4 seconds. The output of the ratemeter was simultane- 
ously recorded, along with relative laser power on the chart recorder. The digital 
recording imit consisted of a gated counter which displayed the pulse count at 1-second 
intervals and a BCD printer which recorded these 1-second count rates. 


Experimental Techniques 

Once the optical components were alined and the detecting electronics were condi- 
tioned for optimum pulse counting, the operation of the equipment was essentially rou- 
tine, requiring only minor fine-tuning adjustments. Prior to each test, the electronic 
equipment and the laser were operated for a 30-minute warmup period. This was neces 
sary for both stabilization of laser power and reduction of spurious counting noise. Dur 
ing this period, the gas cell was purged and filled with nitrogen and the gas temperature 
and pressure were adjusted. 

With the system filled and stabilized, the optical components were adjusted for an 
optimum Raman signal. The rotational Raman spectra were then obtained by scanning 
from the nitrogen anti -Stokes J=20 line at 4839 A upward in wavelength through the 
Stokes J=20 line at 4921 A. The spectrum shown in figure 2 is only a portion of the 
typical scan range. All tests were conducted with room lighting at a minimum to reduce 
the bacl^round count rate. During the tests, all system variables were continuously 
monitored and recorded. Although the monochromator wavelength scale was calibrated 
periodically throv^hout the investigation, no attempt was made to precisely determine 
the Raman line positions and those shown in figure 2 may be in error by ±0. 5 A . 
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TABLE 1. - STANDARD DEVIATION OF 


CALCULATED TEMPERATURES 


[Temperature, 300 K.^ 


Raman lines 
in ratio 

Total Raman intensity of J = 6 
Stokes line, lj_0i counts/sec 


^1 

^2 

50 

100 

300 

500 




Standard deviation of tern- 




perature, a,p, K 


(S)0 

(S)5 

354.0 

250. 3 

144. 5 

111.9 



(S)10 

89.6 

63. 3 

36.6 

28.3 



(S)15 

57. 7 

40.8 

23. 5 

18.2 


■ 

(S)20 

43.3 

30.6 

17.7 

13.7 

(S)l 

(S)4 

566.3 

400. 4 

231. 1 

179. 1 

(S)l 

(S)10 

96. 7 

67. 7 

39.0 

30. 3 

(S)l 

(S)16 

46. 7 

33.0 

19.0 

14.8 

(S)2 

(S)6 

195. 5 

138.2 

79. 8 

61. 8 

(S)2 

(S)12 

52. 7 

37. 3 

21. 5 

16. 7 

(S)2 

(S)18 

38.2 

27. 0 

15. 6 

12. 1 

(S)3 

(S)12 

62. 7 

44. 3 

25.6 

19. 8 

(S)4 

(S)8 

124. 8 

88. 3 

50. 9 

39. 5 

(S)4 

(S)16 

39. 3 

27. 8 

16.0 

12. 4 

(S)5 

(S)20 

44. 3 

31. 2 

18. 1 

14. 0 

(S)6 

(S)16 

42.6 

30. 1 

17. 4 

13. 5 

(S)6 

(S)18 

41, 3 

29. 2 

16. 9 

13. 1 

(S)6 

(S)20 

44. 1 

31. 2 

18.0 

13.9 

(S)7 

(S)18 

46. 1 

32. 6 

18. 8 

14.6 

(S)8 

(S)12 

86. 7 

61. 3 

35. 3 

27. 4 

(S)9 

(S)16 

60. 4 

42. 7 

24. 6 

19. 1 

(S)10 

(S)20 

54. 2 

38. 3 

22. 13 

17. 1 

(AS)0 

(AS) 5 

225. 5 

159. 4 

92. 0 

71. 3 

(AS)0 

(AS) 15 

56. 7 

40. 1 

23. 2 

17. 9 

(AS)2 

(AS) 12 

48. 8 

34. 5 

19. 9 

15.4 

(AS)2 

(AS) 18 

43. 4 

30. 7 

17. 7 

13. 7 

(AS)4 

(AS) 16 

43. 2 

30. 5 

17. 6 

13. 7 

(AS) 6 

(AS) 16 

47. 8 

33. 8 

19. 5 

15. 1 

(AS)6 

(AS) 18 

49. 1 

34. 7 

20. 0 

15. 5 

(AS)7 

(AS) 18 

54. 2 

38. 3 

22. 1 

17. 1 

(AS) 8 

(AS) 12 

90.0 

63. 6 

36.7 

27. 9 

(AS) 10 

(AS)20 

69.2 

48. 9 

28. 2 

21. 9 

(S)0 

(AS) 15 

55.2 

39.0 

22. 5 

17. 5 

(S)2 

(AS) 12 

44. 5 

31. 4 

18. 1 

14. 1 

(S)2 

(AS) 18 

41.6 

29.4 

17. 0 

13, 2 

CS)4 

(AS) 16 

39. 7 

28.0 

16. 2 

12. 5 

(S)6 

(AS) 16 

42. 3 

29. 9 

17. 3 

13.4 

(S)7 

(AS) 18 

48. 1 

34.0 

19.6 

15.2 

(S)8 

(AS) 12 

61. 9 

43. 7 

25.2 

19.6 

(S)10 

(AS)20 

60, 5 

42. 8 

24. 7 

19. 1 

(S)12 

(AS)2 

59. 9 

42. 4 

24. 4 

18. 9 

(S)12 

(AS)8 

170. 1 

120. 3 

69. 4 

53.8 

(S)15 

(AS)0 

59. 8 

42. 3 

24. 4 

18.9 

(S)16 

(AS) 4 

44. 1 

31.2 

18.0 

13. 9 

(S)16 

(AS)6 

50. 5 

35. 7 

20.6 

16,0 

(S)18 

(AS) 2 

40. 3 

28. 5 

16. 4 

12. 7 

(S)18 

(AS)7 

54. 3 

38. 4 

22.2 

17. 2 

(S)20 

(AS) 10 

65.0 

46.0 

26. 5 

20.6 



TABLE II. - OPTIMUM TEMPERATURE 


MEASUREMENT CONDITIONS 





Temperature, 

T, K 





250 



300 



315 




=- M_ K 










*^1 

^2 

M 

Ji 

J2 

M 

Jl 

J2 

M 

(S)2 

(S)16 

227. 1 

(S)2 

(S)18 

270.0 

(S)2 

(S)18 

281. 4 

(S)4 

(S)16 

233.3 

(S)4 

(S)18 

274. 3 

(S)4 

(S)18 

284. 8 

(S)2 

(AS) 14 

233. 5 

(S)2 

(S)16 

274. 3 

(S)2 

(AS) 16 

287.0 

(S)2 

(S)14 

237.8 

(S)2 

(AS) 16 

275. 7 

(S)2 

(S)16 

290. 9 

(S)4 

(AS) 14 

240. 3 

(S)4 

(S)16 

277. 9 

(S)4 

(AS) 16 

291. 1 

(S)2 

(S)18 

240.6 

i(S)4 

(AS) 16 

280. 4 

(S)4 

(S)16 

293. 7 

(S)16 

(AS) 2 

243. 3 

(S)2 

(AS) 14 

281. 4 

(S)2 

(S)20 

293. 9 

(S)2 

(AS) 12 

245.8 

(S)18 

(AS) 2 

285.0 

(S)18 

(S)2 

297. 3 

(S)2 

(AS) 16 

246.0 

(S)4 

(AS) 14 

285.7 

(S)2 

(S)14 

298.2 

(S)3 

(S)16 

246.7 

(S)2 

(S)20 

289.9 

(S)4 

(AS) 14 

301.6 



TABLE m. - STANDARD DEVIATION OF TEMPERATURE CALCULATED 


FROM RATIO OF SUMMED LINE INTENSITIES 


[Stokes -to -Stokes even lines; temperature, 300 K; = 500.] 


Raman line intensities in ratio 

Standard de- 
viation of 

Raman line intensities in ratio 

Standard de- 





viation of 

Numerator 

Denominator 

temperature, 

Numerator 

Denominator 

temperature, 



CTrp, 



CJrji, 



K 



K 

iS + I® 

‘2 + h 

I® 

62.7 

T® 4. I® + I® 

n n ^6 

iS + tS 

h *' ^10 

20.7 



30.6 


ho + h2 

13.6 


I^ 

no 

19. 5 


+iS 

h2 + h4 

10.3 


i^ 

n2 

14.4 


iS +tS 
h4 he 

8.8 


I® 

M4 

12.0 


iS +tS 
he he 

8.3 


I® 

he 

Q 

11.0 


he * *20 

8. 5 


ns 

no 

11.0 

12.0 

yS , tS jS tS 

I2 + I4 + 10 + Iq 

I® +1® 

ho + h2 
h2 + h4 

14.7 
10. 8 






jS jS jS 

12 + 14 + 10 

*8 

34.8 


I® +1® 
h4 + he 

9. 1 


ho 

20 . 5 


he ^ he 

8. 5 


h2 

14.8 


I® +1® 

he + ho 

8.7 


jS 

M4 

ne 

rS 

12,2 

11.2 

4- + yS 

^2 M ^6 

*8 + ho + h2 
ho + h2 * h4 

15, 1 
10.7 


hs 

jS 

no 

11. 2 
12.2 


I® +1® + 1® 
h2 + *14 * ‘le 

I® +1® + 1® 

*14 * *ie + *18 

8. 5 
7. 5 






tS I T® I t® . t® 

‘2 M *6 *8 

iS 

no 

23,3 


iS .s ,s 
*ie ■* *18 * *20 

7.2 


0 






n2 

n4 

ne 

ns 

no 

15.9 
12. 8 
11.6 
11.6 
12.6 

tS ^ yS 4. tS ^ tS 

iS +tS ^tS 
*10 ■* *12 ■* *14 

I® +1® + I® 
*12 * *14 * *ie 

I® + 

*14 ■* *ie *18 

I® +1® + T® 
*16 + *18 *■ *20 

12. 5 
9. 5 
8. 1 
7.6 




T® + 4 - + T® 

*10 ■* *12 ■* *14 ■* *16 

tS . tS tS jS 

*12 *14 *16 ■* *18 

I® +1® +1® + T® 
*14 * *16 ^ *18 ^ *20 



jS . tS 

*2 * M 

i^+i® 

h + h 

34. 5 

if+if+if+if 

9. 4 


h * ho 

19.9 


7 . 7 


ho * h2 

13.6 


6. 9 


jS jS 

h 2 + h 4 

10,5 





h4 + he 

8.9 





tS ,s 
he *18 

8.3 





h8 ho 

8.4 
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